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Detection of Geometric Structures by an Optimally 
Subsampled Shearlet System in Noisy Digital Images 

Philipp Petersen 


Abstract 

We provide a statistical analysis of the ability of digitized continuous shearlet systems 
to detect objects embedded in white noise. We analyze the possibility to subsample the 
shearlet system and obtain a system of significantly reduced cardinality that can still yield 
statistically optimal detection results. 


1 Introduction 

Detection and classification of singularities in images is an important task in imaging appli¬ 
cations. It has been shown in that the continuous shearlet transform can be applied 

to identify location and position of jump singularities of an image. These works have been 
extended and improved in [9] to include compactly supported shearlets, which have optimal 
spatial localization and can therefore handle singularities better. 

All of the above results use a continuous shearlet system i.e. a system constructed from a 
generator p € L 2 ( R 2 ) together with elements 

ipa,sA x ) = - *)), a € (0,1], s € [-1,1], t € M 2 , (1) 

where 

a “=(o D “ d s * = (J 0 - 

For applications we need to sample the parameters a,s,t and also restrict the functions 
ip a ,s,t to a digital lattice or an n x n pixel image. 

In this work we will examine this digitization from a detection point of view, in the sense, 
that we aim to describe a suitable digital shearlet system to optimally detect geometrical 
structures embedded in noise. This is inspired by the work of Arias-Castro, Donoho and Huo 
[2], which describes how to construct statistical tests to test for some geometric structures 
such as lines, intervals or rectangles. There it is assumed that £ is a discrete signal of length 
n E N with certain characteristics such as jumps and a white noise signal z(i), where i is an 
element of some discrete index set I n and z(i) ~ A/"(0,1). In the one dimensional situation we 
choose I n := {1,..., n}. Similarly in higher dimensions we let I n := {i = (ii, 12 , ..., ik) ■ 1 < 
il < n, l = 1,... ,k}- The measured signal is now given by 

x(i) = £(i) + z(i), i€l n - (2) 

In |2] this model was used to detect discrete lines, rectangles, circles and intervals embedded 
in white noise. We will continue on this path and examine the optimal detection of structures 
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connected to discontinuities and continuous shearlet systems. Our main result, Theorem 13.61 
states, that careful subsampling of a continuous shearlet system does not influence its abilities 
to detect certain geometric properties. This subsampling yields a reduction of the size of the 
system from 0(n 4 ) to 0{n 2 log(n)) for an n x n pixel image. The manuscript is structured 
as follows. We will describe the methodology of using statistical tests for the detection of 
geometrical structures in Subsection 12.11 Then we will proceed to construct a test to find 
jump sdiscontinuities in Id signals in Subsection 12.21 From these methods we will extract a 
general framework in Subsection 13.11 This will then be applied to the continuous shearlet 
system in Subsection 13.21 

2 Detection of Jump Singularities 

2.1 Problem Setup 

We will begin with the detection of jump discontinuities. From the point of view of detection 
this geometrical structure is closely related to that of intervals and therefore we will first 
briefly recall the detection of intervals. Using the model Q with £ = A n £r a w (i), where 
£[a,bl(*) = Vfr — a X[a,b], In = 1 , ■ ■ ■ n, the authors of [2] showed, that a test can reliably the 
presence of an interval in an asymptotically stable way, that will be introduced below, from 
the signal 


x(i) = A n ^ b] (i) + z(i), 1 <i<n, (3) 

if A n lies above a critical threshold. On the other hand if A n lies below this threshold, no 
test can reliably detect the interval. The phrase of asymptotically reliable detection needs 
some mathematical rigor and hence, we repeat the definitions of asymptotically powerful and 
asymptotically powerless tests. 

Definition 2.1. JBj In a sequence of testing problems (Uo.n) vs. we sa V that a sequence 

of tests (T n ) is asymptotically powerful if 

Ph 0 ,„ {T n rejects H 0 } + P Hl n {T n accepts H 0 } 0, n ->■ oo 

and the sequence is asymptotically powerless if 

Ph 0 ,„ {T n rejects H 0 } + n {T n accepts H 0 } ->• 1, n ->• oo. 

2.2 Detecting ID Jump Singularities 

We aim to construct a series of tests that are asymptotically powerful, to detect jump singu¬ 
larities in noisy data. This means, that we obtain data of the form 

x(i) = £j(i) + z(i), 1 <i<n, (4) 

where £j exhibits some sort of jump, and z(i) ~ A£(0,1). In the spirit of [2] one must first 
pose the question under which conditions any test can reliably detect the presence of a jump 
in the signal. 

First of all we need to determine what a jump singularity should be. At this point we may 
want to consider the prototype of a jump singularity, given by objects of the form if(i) = 
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X[l,s](®)>® £ I n which model a jump at the position s from 1 to 0. These types of jumps 
are not very challenging, in fact the detection of these is already covered by the possibility 
to detect intervals as in ([3]). Another approach could be to model jumps simply as a point 
where the function value differs strongly from that of the neighboring points. Taking pairwise 
differences leads towards the problem of finding one elevated mean within white noise, i.e. the 
needle in a haystack problem, El HZl- In order to be reliably detectable this would require the 
jumps to be quite high and in fact higher than we might want to. In fact, a jump could also 
be detectable, if it is small but in a neighborhood apart from the jump the function values of 
the signal change very little. In other words, if before and after the jump we have an almost 
constant signal, then we would also expect to be able to detect smaller jumps. Indeed, a 
reasonable model for a jump singularity is that a function has a large Haar wavelet coefficient. 
This will now be our model. 

Let us examine the continuous Haar wavelet transform on 1 < i < n. Let for 1 < a < n/2, 
1 < t < n, 

i>a,t(i) = -7=(X[0,a](i-t + a)-X[a,2a](' i - t + a )), 1 < * < (5) 

\J la 

where the expression i — t + a is to be understood modulo n. We denote by 

W := ^ < t < n,l < a < n/2] 

the continuous Haar wavelet system. The cardinality of the continuous Haar wavelet system 
is n 2 / 2. We choose the continuous Haar wavelet system as a model for jumps in the signal 
i.e. a signal has a jump of size A at t, if | (£,i/j a ,t) \ = A for some 1 < a < n/2. 

Of course our jump singularity model is not the only possibility. One could certainly also 
work with unbalanced jumps, which means that we only need a smooth part with small 
variation on one side of the jump. Should one take the continuous Haar wavelet transform of 
white noise, i.e. look at the values a simple calculation yields that we obtain n 2 /2 

values that are normal distributed. If the values {i/> a) t,x) where to be independent, the needle 
in a haystack problem suggest, that detection of a jump is only reliably possible if the jump size 
is larger than \J2 log(n 2 ). However, the probability variables x ) are not independent and 
we will show that the optimal threshold for reliable detection is sup a t \ {ip a ,u x ) I > \/2 log(n). 
This means, that if for some rj > 0, sup at | ( 4> a ,t , x ) | < (1 — r/)y21og(n) there does not exist 
an asymptotically powerful test, but if sup at | {if at t,x) | > (1 + r})y/2\og(n) we can construct 
an asymptotically powerful test. 

The key to show that these lower thresholds can be achieved is to subsample the scaling 
parameter of the continuous Haar wavelet system. For e > 0, we choose values a, = 2 ek ,k = 
1,... [^(log(n) — 1)] and define 

W e := | il>a k ,t,k = 1,... |"^(log(n) - 1)],1 < i < n| . 

We use the metric: ^(^ 1 ,^ 2 ) := (1 — (£ 1 , £ 2 ))^ to describe the difference between the elements 
of and the continuous Haar wavelet system W. The following Lemma shows, that the 
subsampled system is in some sense dense with respect to the metric 5. 

Lemma 2.2. For e > 0 the subsampled Haar wavelet system W e is an e-net for the continuous 
Haar wavelet system W, i.e. 

max min 6(if ak , s > Vw) < e- 

a,t k,s 
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Proof. For given a < nj 2, there exists > 0 such that either a& < a < 2 e / 2 afc or 2 e / 2 a£ < 
a < ak -|-i • Then 

i ^aktt) — 1 (Vki,t) 'lpai e ,t) > 

where 


(VVi,i) lfa k ,t) — (lfa k ,ti l(a k ,t) T 'f^a,tt'^fa k ,t) — 1 T (lfa k ,t VL,i j f^a k ,t) ■ 

Furthermore 

(lfa k ,t "0a, tj lfa k) t) — IIV’afc,t \ \ \ | ^fa k ,t \ \ "Pi HV’a/ct Vb,t||- 

Assume a& < a < 2 e// 2 afc, then applying the definition ([5]) we see that — ?/a,t || 2 can be 

estimated by 277 -la — ad + aJ- 7 =-J=| 2 . Now we have 

•^ a \/2 a v 

|a - a fe | < a fc ( 2 e/2 - 1 ) < e/ 2 . 

This yields that 

|1 - (lf a ,t,lfa k ,t) | < e 2 

and hence 5(if a ,t,ipa k ,t) < e. 

The case of 2 e / 2 afc < a < a^+i follows analogously. □ 

We denote by 


W(x) = SUp (| (ljja,t,x) \,lfa,t € W)} 
the maximum value of the Haar coefficients of a signal x. Similarly 

W e (x) = sup (| (ifa t t,X ) I ,lfa,t € W e )} 

denotes the maximum Haar coefficient of x with respect to the subsampled system W e . Using 
the e— net of Lemma O we aim to now construct a reliable test. Indeed we will show, that 
if the jump size is larger than y / 2(l + 77 ) log(n) for some positive 7 / the generalized likelihood 
ratio test (GLRT), which rejects the null hypothesis if 

W(x) > t n , 

is asymptotically powerful, if we choose t n /y/2\og(n) —>• 1. To prove this, we employ the fact 
that the subsampled Haar wavelet system yields an e-net. Furthermore, we make use of the 
GLRT with the reduced dictionary, i.e. we reject the null hypothesis if 

W e (x) > t n , (6) 

where t n j yj2 log(n) —>• 1. The last ingredient for our desired result is the following result 
describing extreme values of Gaussian processes. 

Lemma 2.3. J2j/ Let ui,... ,L 0 m be (possibly dependent) A^(0,cr/) variables with all <jj < er. 
Then, 


P(max{wi,..., Lo m } > v / 21og(m)o- < —=====. 

yf 47T log(m) 
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Now we will prove that asymptotically the reduced GLRT and the full GLRT behave in the 
same way. For this, let 

d e (x) : = min |W(x) - (ifa k ,t,x) I, 

4>a k ,t&We 

denote the effective distance of VV to W e . 

Lemma 2.4. For each r/ > 0, 

p H 0 , n > V81og(n)e} < 1 

1 J Vlog(n) 

Proof. First of all we assume that the maximum of the full GLRT is attained at ip a ,t- By 
Lemma 12.21 we have that there exists f € W e such that || if a ,t ~ H ; a k t Ik < e - Then 

<5(x) < \{^ ak f-^a,t,x)\ = | ^ (V’ofc.tW -^a,t(*))*(*)l- 

l<i<n 


Since we assume that iLo,n holds, we have that x(i) = z(i) ~ IV(0, 1) for all 1 < i < n. 
Now “ ^a,*(*))*(*) ~ Af(0, ||^ afci f(i) - Vvll 2 ) = -^(0,e 2 ) with e < e. To estimate 

the probability of S € (x) exceeding 4ylog (n)e we have to take the maximum of all probability 
variables X^(Va fe f(*)~~ ^a,t(i))z(i) over all a, t and corresponding a*,, t. These are certainly less 
than n 4 variables which are all normal distributed with variance less than e 2 . Hence Lemma 
12.31 implies that 


p Ho,n > V21og(n 4 )e} < , 

1 > yj 47T log (re 4 


1 

< 

V l °s(n) 


□ 

Now we are in the position to prove that the GLRT is asymptotically powerful, should the 
jump size exceed y/2(l + rj) log(n). 

Theorem 2.5. For each rj > 0, 

Ph 0 |w(x) > v/2(l+^)Iog(n)| ->• 0, for n oo. 

Proof. Let us first prove that, for e > 0, 

Ph 0 |>V e (x) > y/2{l + rf) log(n)| -3> 0, for n —> oo. 

In fact jfW e < ^n(log(n) — 1). Since under H o, (if ak ,t,x) ~ jV(0,1) for all if ak ,t € W e Lemma 
12.31 implies that 

Ph 0 |w e (®) > y // 21og(^n(log(n) - 1)) j = P Ho |w e (x) > y/2(\og(n) + log(log(n) - 1) - log(e)) 

1 

y/47rlog(±n(log(n) - 1)) 
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Furthermore we obtain that for any r/ > 0, there exists M v > 0 such that for all m > M v we 
have ?/log(?n) > log(log(m) — 1) — log(e). Hence 

Ph 0 |w e (x') > y/2{\ + rj) log(m)| < P Ho |>V e (x) > ^(log^) + log(log(n) - 1 ) - log(e))| 

< — = —»■ 0, for n —»■ oo. 

\/ 47 r log(n) 

The result for W(x) follows by choosing e such that ^/2(1 + 77 ) — %/ 8 e < ^/2(1 + 77 / 2 ), then 
Lemma 12.41 and the observation that 

w e (z) < vw(s) < w e (x) + <y e (0, 


imply the result. □ 

Of course the question arises weather the bound of Theorem l2.5l is tight, i.e. if ^/2(1 + 77 ) log(n) 
is the minimum jump size that can be detected. The following Lemma states that this is indeed 
the case, i.e. that a lower jump size will not be detectable with an asymptotically powerful 
test. 

Lemma 2.6. Given a sequence of hypotheses 

H ln '■ = £/(*) + z(i), 


and rj > 0 with 

W e {£,j) = \/2(T - rj) log(n). 

Then, no test is asymptotically powerful for testing iLo,n against H^ n . 

Proof. Let Ho <n be the hypothesis that x{i) ~ jV"(0,1). Furthermore let H[ n be the hypothesis 
that for any 0 < t < n/2 — 1 

x(i) = A n if 1>2 t + z(i),Zi ~ AA(0,1), 0 < i < n - 1. 

Consider the orthonormal Haar wavelet basis on 1 < i < n consisting of 

,2t ■ t G [0, n/2]} U 

Indeed, since this yields an orthonormal basis for the signals on I n we can apply a basis 
transformation and represent x with respect to the Haar wavelet basis to obtain x. Now H[ n 
becomes 




x(i) = A n 6 t (i ) + z(i), 0 < i < n. 

The needle in a haystack problem, suggests, that if A n = \/2(l—r/log(n)) for i] > 0, then 
H() n and H[ n merge asymptotically, i.e. there does not exist an asymptotically powerful test 
to distinguish between the two hypotheses. Since distinguishing H[ from iLo,n should be 
easier than distinguishing Hi <n from iLo,n this proves the claim. □ 
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3 General Geometrical Structures 


3.1 General Framework 


In this section we will extract the main ideas of the path we followed in the previous subsection 
on Haar wavelets to give a theory for arbitrary systems of functions. This can then be employed 
in any dimension and especially in Section [3.21 to two dimensional functions. We saw in the 
Haar wavelet case, that subsampling the functions dyadically does not destroy the detectability 
of jumps by that system. Indeed, we will show, that whenever we have set of functions on I n : 
(gj)jeJu with some index set J n such that for every e > 0 there exist subsystems, that yield 
an e-net the reduced GLRT is asymptotically powerful if the signal has a detectable feature: 

Definition 3.1. Given a function system, ( gj)j£j n with a subsystem , C J n on I n . 

If for each l E J n we have that minjgje \\gj — gi\\ < e, we say that ( gj)j£j * constitutes an e-net 
for {gj) je j n . 

Definition 3.2. For each n € N let ( gj)j£j n be a system of functions on I n with = 

\gj{i)\ 2 = 1- Furthermore assume that for every e there exists a subsystem ({gj)jej%) 
hat constitutes an e-net for (gj)jeJ„- Assume that there exists e > 0 such that for every u > 0 
there exists a constant C u with 

\Jn\< C w n e+UJ , for all n € N. 

A sequence of signals has a detectable feature if sup ?G j n | n \gj ) | > (1 + rj)y/2log(n e ) 

for some r] > 0 and n E N. 

The jump singularities from the preceding subsection are detectable features of the con¬ 
tinuous Haar wavelet system W. We will study another examples in the sequel which are 
detectable features of continuous shearlet systems. Moving on with the general set-up and a 
set of functions ( gj)jej n , we define 

Q n {x) := max | (gj,x) |, 
ieJn 

and for the subsystems we write 


Q e n {x) := max | (gj , x) |. 
jeJ' 

Q n {x) and Gn(x) denote the maximum response from a signal, when tested with Q n or Qf 
respectively. Using these we can again build a full and a reduced GLRT by checking, whether 
G n {x) or Gn(x) exceed a critical threshold. Furthermore, we aim to obtain an analog of Lemma 
12.41 which would tell us that asymptotically we can replace the full GLRT with the reduced 
GLRT. With 


5 e (x) := min | G(x) - (gj ,: 


we obtain 


Lemma 3.3. For each n E N let ( gj)j£j n be a set of functions on I n . Assume that for every 
e > 0 we have subsystems ((gj)j£j^) that yield an e-net for (gj)j^j n - Then, we have 


p H 0 ,„ {<5 e (x) > V81og(n)e} < 


1 

\/ lo g( n ) 
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Proof. The proof is the same as that of Lemma 12.41 


□ 


Now we can show, that every detectable feature can, as the name already predicts, be 
detected. Furthermore it is possible to construct a test using only a reduced GLRT. 

Theorem 3.4. For every n € N let ( gj)j£j n be a set of normalized functions with subsystems 
as in Definition 1,9.111 that have cardinality 0(n e+u: ) for e > 0 and every uj > 0. Let £f n ’ be a 
sequence of signals that have a detectable feature and 

x n (i) = + z(i), where z(i) ~ Af( 0, 1 ) and i € I n . 

Then there exists e > 0 such that for all rj > 0 

P | Q e {z) > (i/ 2 (l + r 7 ) log(n e ))| ->• 0 , for n —00 

and there exists 77 > 0 such that 

P j<5 e (x < ' n / > ( 1 / 2(1 + t?) log(n e ))| —> 1 , /or n —>• 00 . 

Proof. The first part of the proof uses similar arguments as the proof of Theorem 12.51 For 
e,rj > 0, there exists 00 < 77 and C u such that ffQ e < C u n e+U1 . Since the gj are normalized we 
obtain that ( gj,z ) ~ jV(0,1). Lemma [2.31 implies 

P {w e (z) > V / 21og(C 0 J n e +-)} = P [W e (z) > ^2(1 + w) log(ra c ) - log^))} 

< 1 
~~ y^vr log(C^n e+a; ) 

Furthermore we obtain that there exists M r] > 0 such that for n > A4 r] we have 77 log(n e ) > 
u log(?r e ) — log(C^). Hence 

P | W £ (z ) > 1 / 2(1 + T 7 )log(n e )| < P jw e (z) > i/2(l + w) log( 77 e ) - log(C(/)} 

1 

< — . — 0 . 

y 47 T log(?z e ) 

For the second part we first prove that there exists some 77 > 0 such that 

P |<5(x (n) ) > (\/ 2 (l + 77 ) log (n))j ->• 1 , for n — > 00 . (7) 

By the assumption that x has a detectable feature we obtain that there exists some 7 such 
that Q( x( n )) Af(y/2(1 + 7) log(n e ), 1). Choosing 77 < 7 less than 7 yields that (jTj) is equal 

to 

P |A/(v / 2(l + 7 ) log(n e ), 1 ) > i/ 2 (l + 77 ) log(n)| = P |jV(0,1) > -fiy/logn^ < 0 , 


where the last equation follows by Mills ratio and (5 = -/2(1 + 77 ) — ^2(1 + 7). To obtain 
that 

P{e e (x (n) ) > (i/ 2 (l + 77 ) log(n))| 




















we observe that 



Let g-j £ Q be a function where , then 



The rest follows by Lemma 13.31 


□ 


Let us reflect shortly on what we proved so far. If a signal has a detectable structure, 
i.e., a structure that can be detected by a set of functions, we can subsample this set 
to obtain a cheaper test. We saw analyzed how this can be applied for Haar wavelets. 
Jump singularities are detectable by the continuous Haar wavelet system as we observed 
in Subsection 12.21 Furthermore the continuous Haar wavelet system can be subsampled 
to yield e-nets as shown in Lemma 12.21 This smaller set allows for the construction of a 
fast test by computing all the scalar products of the subsampled Haar wavelet coefficients 
with the signal. In fact this can be done in 0(nlog(n)), which we see by observing, that 
(V’afc.t) = ip a ,o * x(t) = ifft(fft(V’a,o)fft( a: ))) where fft and ifft are the discrete fast Fourier 
transform and its inverse. Since we have only 0(log(?r)) scaling parameters, this yields the 
claimed computational complexity of 0(nlog(n)). Since this machinery seems quite powerful 
we can now advance to higher dimensional structures namely those described by continuous 
shear let systems. 

3.2 Detecting 2D Jump Singularities 

From the work of [9] is is clear, that continuous shearlet systems are well versed to detect 
geometric structures like discontinuities along smooth curves. By virtue of Theorem 13.41 we 
need to show, that jump singularities are a detectable feature of a continuous shearlet system 
and that subsampling the shearlet system yields a subsystem, that constitutes an e-net. 

Let us start by defining the continuous digital shearlet transform onanx n pixel image. As 
the digital realm should be a discretization of the continuum realm L 2 ([0, l] 2 ) this should be 
reflected in the construction. Let ^ be a continuous shearlet supported in [0, l] 2 . We define 
the continuous digital shearlet system as digitization of the continuous system by defining for 
an n x n pixel image the continuous digital shearlet system CT>ST~L as 


CVSU := 



(h,t 2 ), 1 < U < nj , 


where 
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In the equation above Q q (r) denotes the cube of sidelength q and center r. As a second 
step we define subsystems with significantly lower complexity that still yield e—nets. With 
[~£~| = (fti], [£ 2 !)) we define for rpu,v > 0, 


CDSTL 5 ^'" := ak = 2&k ' s = h = ua k r i, t 2 = vjaf/nr 2 , 


k = 0 , 


-(log(n) - 1) 


,l = - 

y/n 


y/n 

_ yjQ'kLd _ 

, . . . , 

_ yj akU_ 

n 

0,..., - 

, , r 2 = 

= 0.... 

Vn \ 

a k u 

y/a k i / J 


Observe that the cardinality of C'DS'H S ’ U},U is of order n 2 (log n ), which compared to the colossal 
0(n 4 ) of the CDSPL can be seen as a significant reduction. In order to apply Theorem 13.41 we 
still need to see that the new system constitutes an e-net. 


Lemma 3.5. Let ip be a Lipschitz continuous shearlet with Lipschitz constant C^. Then for 
S < 1. C'DSTL 5,U} ’ 1 ' is an e-net for CVSTL with e = 25/3 + C^35 + oj + 3 C^v. 

Proof. We compute 


where we need to choose k, s and f in a suitable way. Without loss of generality, we can 
assume t = 0, since the above expression is invariant under translations. Furthermore, we 
choose 


a k = argmin | \ai/] — a|, s= argmin |s — fs]|. 

k §=—Tj- :ujy/a^n: ^ 

Assume for simplicity that a k < a. Then we obtain the error estimates 

I \a k ~\ — a\ < a k (2 d — 1), \s - \s\\ < ujy/a/fh/2. 

Moreover,choose t = (t\,t 2 ) such that 

~ v ~ v _ 

|£i| < tj a k , |£ 2 | < -yja k n. 

Then we obtain 


^a,s,0 ,s,t II 
= ka, S , 0 (*) - V’afc.S.tM 


i£.In 


E 

E 


5 

n4 


5 

n4 


/ 

JxGQ i 

2 n 

/ 

Jx£Q i 


(i/n) 


a iip{Aa l Ss_ 1 x) - \a k \^ip{A l f h] S^{x-i/n))dx 

n n —— 


-i c-i 


a Zip(Aa 1 S ± 1 x) - \a k ] S - A ^ S^i/n)dx 


kl 

n n 


-1 o—l, 


t-1 o-l- 


(i/n) 


n n 
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This implies that I 2 < 2l 5 C^\x-^x\+ {s ^ ) J° \x 2 \). Using that = 2 5 and ( %[^ < 

uj/2 we obtain that I 2 < C^,2^ ((2 s — 1) + w/ 2 )||x ||2 and since ||a :||2 < \/2 we get I 2 < 
C^ S+ ^((2 S - 1) +w/2) < C^S6 + uj. 
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With a similar approach we can estimate I 3 by 


2 



t/n = 2 4 5 C^ 


< 2 \ S C, 




y/n\ak\ 


ak v ytTkls] V y/ak V 

fofcl 2 Vn\a k ] 2 ^\a k ] 2 


< 3 C^v. 


In conclusion we obtain that 


Ii +12 + -I 3 ^ 25/3 + C^35 + w + 3Cipu 


which yields the claim. □ 

Next we can combine Theorem EH and Lemma 13.51 to obtain, that subsampled shearlet 
systems can detect all detectable features of unsubsampled shearlet systems. 

Theorem 3.6. Let if: be a Lipschitz continuous shearlet and CDSTL the corresponding contin¬ 
uous digital shearlet system. Let £( n ) be a sequence of signals that have a detectable feature. 
Then, there exists 5, lo,v > 0 such that the GLRT with CDSTL s,u1 ’ 1 ' is asymptotically powerful. 


3.3 Open Questions 

What can be detected by CDSTL ? 

First of all it is not yet clear, what type of geometrical structures can be detected with the 
CD STL. Certainly, elements of CDSTL are detectable with the GLRT and a reduced shearlet 
system CDSTL S,UJ,U , but one would certainly be interested in objects or structures that appear 
naturally in images. Since the publications EJEliSj show how shearlets can detect jump 
singularities, an extension to these structures seems very likely. 

Is the subsampling optimal? 

We see that the cardinality of the subsampled shearlet systems CDSTL 5,UJ ’ 1/ is 0(n 2 log(n)). 
Given that the image has n 2 pixels, it seems very likely that, at least up to the log factor, this 
is the minimal cardinality for which an e-net can be constructed. A proof for this optimality 
is still missing. A different approach would be to allow further subsampling on the cost of an 
increasing signal strength. If this is possible, a precise characterization of the trade-off as in 
[5] is desirable. 

Extensions 

Similar to the question of what can be detected, we can also adjust the process of testing. 
For instance if we assume we have an image composed out of only a few elements of CDSTL, 
we can develop a procedure, where we find some promising elements in CDSTL S ' u,l/ and try 
various combinations of them. A similar method to extract curves from noisy images based 
on extensions of promising parts of curves has been proposed in |TJ. Since shearlets yield 
optimally sparse approximations of cartoon-like images [8], this might ultimately give rise to 
a method to optimally detect cartoon-like images embedded in noise. 
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